The characteristic blue spectra of accretion disks in quasars 
as uncovered in the infrared 
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Quasars are thought to be powered by supermassive black 
holes accreting surrounding gas^"^. Central to this picture is 
a putative accretion disk which is believed to be the source 
of the majority of the radiative output^"^. It is well known, 
however, that the most extensively studied disk model^ — 
an optically thick disk which is heated locally by the dissi- 
pation of gravitational binding energy — is apparently con- 
. tradicted by observations in a few major respects*'^. In 
particular, the model predicts a specific blue spectral shape 
asymptotically from the visible to the near- infrared^*, but 
this is not generally seen in the visible wavelength region 
where the disk spectrum is observable'"'^. A crucial diffi- 
culty was that, toward the infrared, the disk spectrum starts 
to be hidden under strong hot dust emission from much 
larger but hitherto unresolved scales, and thus has essen- 
tially been impossible to observe. Here we report observa- 
tions of polarized fight interior to the dust-emiting region 
that enable us to uncover this near-infrared disk spectrum 
in several quasars. The revealed spectra show that the near- 
infrared disk spectrum is indeed as blue as predicted. This 
indicates that, at least for the outer near-infrared-emitting 
radii, the standard picture of the locally heated disk is ap- 
proximately correct. 

A success of the most extensively studied disk model is 
that it gives the radiative output peak approximately correctly in 
the ultraviolet (UV; ~0. 01-0.4 //m) wavelengths for the case of 
a supermassive black hole. This is observed for a generic spec- 
tral energy distribution'"* of quasars, the most luminous exam- 
ple of active galactic nuclei (AGN). However, it has long been 
known that the model apparently shows a few major contradic- 
tions with observations^'^. One of the disagreements, and per- 
haps the most easily comprehensible one, is the spectral shape 
of the radiation. From the basic hypothesis of the model (an op- 
tically thick disk heated locally), the effective disk temperature 
T is fixed as a function of radius r as T oc r^-^/'^ over a broad 
range of radii. This leads to a well-known blue spectral shape 
limit, Fi, ex ^z+^^Z'^, being asymptotically reached at long wave- 
lengths from the visible (also called the optical; ^0.4-1 ^m) to 



the near-infrared (^1-2 /im) for AGN disks. In contrast, many 
studies have shown that the general AGN spectral shape ob- 
served at optical/UV wavelengths is much redder, spectral slope 
a being from —0.2 to —1 (Fi, oc j^"), and never as blue as this 
spectral shape^"'^. 

The predicted blue shape limit is strictly true in the sim- 
plest assumption of local black-body emission. In more so- 
phisticated disk atmosphere models'^, the spectrum generally 
becomes slightly redder at optical wavelengths, owing to vari- 
ous opacity and non-LTE effects, but discrepancies between the 
model and observed spectra still remain'^. However, the red- 
der model slopes at optical wavelengths form a wider concave 
spectrum that shifts the bluer limit above to longer wavelengths, 
into the near-IR. The observed spectra certainly appear to be- 
come bluer from the short UV wavelengths to the optical. The 
crucial observational difficulty here has been that the disk spec- 
trum starts to be hidden under the hot dust thermal emission 
which begins at wavelengths greater than ^^1 /im, a limit set by 
the sublimation temperature of dust grains (-^1500 K). These 
dust grains exist at larger spatial scales, in a configuration of- 
ten thought to have a torus-like geometry but which is generally 
not yet spatially resolvable. Therefore, it has been virtually im- 
possible to observe the underlying near-IR disk spectrum'^. We 
note that, contrary to early spectral fitting studies^**'^''^, the 
IR component is no longer thought to be non-thermal, and thus 
cannot be extrapolated to underlie the optical spectrum — an 
extrapolated non-thermal spectrum had effectively made the in- 
ferred disk spectrum on top of the non-thermal spectrum bluer 

We argue here that this buried part of the disk spectrum 
can be revealed by observing the near-IR polarized light. Op- 
tical continua of many directly- visible AGNs called Type Is 
(Seyfert 1 galaxies and quasars) are known to be linearly po- 
larized at a level <1%. The polarization position angle (PA) in 
these Type I cases is mostly parallel to the rotation axis of the 
putative accretion disk, where the axis can be probed by the lin- 
ear jet-like structure of radio emission'^"-". (This is in contrast 
to the cases in hidden AGNs or Type 2s, which show high polar- 
ization at perpendicular R4s and strong broad lines in polarized 
light, and which are not the subject of the present paper.) This 
polarization in Type 1 s is interpreted as an indication of an equa- 
torial scattering region, which is optically thin and surrounds the 
disk. In many Seyfert 1 galaxies, broad emission lines are polar- 
ized at a much lower level than the continuum polarization and 
at different PAs'^'^*', indicating that the scatterers reside roughly 
at the same spatial scales as the broad-line-emitting clouds. 

At least in several quasars, the emission line polarization 
even vanishes — the optical polarized light spectrum shows no 



or very little emission line flux^''^^. This is very likely to indi- 
cate that the scatterers reside interior to the broad-line clouds. 
In these Type 1 cases, the scatterers are thought to be elec- 
trons and not dust grains, since the scattering region is interior 
to the dust sublimation radius. Because electron scattering is 
wavelength independent, the polarized light therefore produces 
a copy of the spectrum originating interior to the scattering re- 
gion. (We note that this electron scattering is conceptually dif- 
ferent from the one in previous works^^-^"^, which was assumed 
to be intrinsic to the accretion disk atmosphere and gave rise to 
the prediction of perpendicular PAs.) In the optical polarized 
light from these quasars, which excludes the emission from the 
broad-line region, we actually found a hydrogen Balmer-edge 
feature in absorption for the first time, which we believe origi- 
nates in the disk and is a very specific indication of the thermal 
and optically-thick nature of the emission^'"^^. 



We can then use the same polarized light, but in the near- 
IR, to reveal the hidden spectrum of the disk, by stripping off the 
dust radiation from the torus exterior to the scattering region. 
In previous work^^, we suggested that this method appears to 
work in at least one quasar that has polarization data with high 
signal-to-noise-ratio at two wavelength bands in the near-IR. 
However, crucial information needed at the time was whether 
the near-IR polarized light behavior is consistent and systematic 
in different objects. If it is, this would critically argue against, 
for example, a possible secondary polarization component re- 
lated to dust grains newly arising in the near-IR. Therefore we 
undertook the near-IR polarimetry of five other quasars. The 
targets were selected to be polarized in optical continua but es- 
sentially not in emission lines, to ensure the scattering to be 
interior to the broad-line region. These properties were either 
already known^^ or were determined in our optical polarimetric 
survey and follow-up spectropolarunetry. 
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Figure 1 : Overlay of the polarized- and total-Ught spectra ob- 
served in six different quasars. We plot scaled vFi, data: Q0144- 
3938 (redshift 2=0.244), green; 3C95 (z =0.616), blue; CTS A09.36 
(2=0.310), light blue; 4C 09.72 (2=0.433), red; PKS 2310-322 
(2=0.337), light green. Plotted in purple are the data for Ton 202 
(2=0.366) from a previous paper . Total-light spectra, shown as bold 
traces in the optical and as squares in the near-infrared, are normal- 
ized at 1/im in the rest frame, by interpolation (except for 3C95 for 
which we normalized by z^F^ observed at 1.3 fim in the rest frame). 
Polarized-light spectra, shown as light points in the optical and as bold 
points in the near-infrared both (vertical error bars, 1-cr), are sepa- 
rately normalized, also at l/im, by fitting a power-law to the near-IR 
polarized-light spectra. For both total-light and polarized-light data, 
horizontal bar lengths indcate bandwidth. The normalized polarized- 
light spectra are arbitrarily shifted downwards by a factor of three rel- 
ative to the normalized total-light spectra, for clarity. The total-light 
spectra in vFi, turns up at around or slightly longward of 1 /^m. In 
contrast, the polarized-light spectra in vFi, all consistently and sys- 
tematically decrease towards long wavelengths, showing a blue shape 
of approximately power-law form. 



In Fig.l we show the spectra of linearly polarized light 
measured in the near-IR broad-band imaging polarimetry, and 
optical spectropolarimetry, of the six quasars in total (includ- 
ing the one studied in the previous work). Details of the mea- 
surements, as well as the procedures for removing instrumental 
polarizations, can be found in the Supplementary Information. 
Generally, polarization degrees observed for these quasars are 
~I% at ~0.5 /xm, gradually decreasing to ^0.5% at ^2 /im, 
and PAs are essentially constant over the near-UV/optical/near- 
IR wavelengths for a given object. A significant result here is 
that all the objects behave in a similar and systematic way show- 
ing blue polarized-light spectra. While the total-light spectra in 
vFi, turn up at around 1 /im toward longer wavelengths due to 
the onset of dust emission, all the polarized light spectra, which 
eliminate dust, display a rapid decrease in i/Fiy, with a shape of 
approximately power-law form. 

The measurement of the spectral index a in F^, (oc t^") 
for each object is shown in Fig. 2. The measured slopes are con- 
sistent with each other within their errors, and the individual 
slopes as well as their average clearly point to a shape much 
bluer than those observed in the UV/optical. Astonishingly, they 



are all consistent with the F. 



y+l/3 



shape. The weighted 



meanof the measured slopes is a = +0.44±0.11. Although the 
sample size is small, there does not appear to be any luminos- 
ity dependence, as seen in Fig.2, and we did not find dependen- 
cies on black hole masses Mbh or Eddington ratios L/L-Edd de- 
rived from the width of the H/3 Balmer line. This is expected if 
the near-IR spectrum is in the long-wavelength limit of the disk 
model, which is independent of parameters such as black hole 
mass or Eddington ratio. In this case, by regarding each mea- 
surement as a measurement of the same quantity, the weighted 
mean over these measurements becomes physically meaningful. 
We note that, if we formally convert the mean slope to the ra- 
dial temperature distribution for the case of an optically-thick 
disk, we obtain T oc 7'-0-78±o.03^ consistent with the predicted 
dependence T ex r^"^/"*. 

The systematic behavior of the near-IR polarized light, 
as well as the constancy of the PAs over all wavelengths. 
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Figure 2: Spectral index of polarized light spectra. We plot a (in 

Fi, oc v") against vL,j for total light at 0.51 ^m. The index was 
measured by a power-law fit for each near-IR polarized-light spectrum 
(note the different wavelength range covered depending on the redshift) 
and is shown with 1-a error bars. A weighted mean of the spectral 
index measurements is shown dashed; the shaded area represents its 
deduced l-a uncertainty. The mean or median slopes of the UV/optical 

9—1 ^ 

total-light spectra derived in various studies are also shown. 



strongly argues against there being any secondary polarization 
contamination. We might worry that the polarized light would 
be affected if the corresponding spatial scale of the disk emis- 
sion at long wavelengths were to become large and finite in 
comparison with the size of the scattering region. However, this 
seems unlikely, since the half-light radius of the disk (within 
which half of the total light is emitted) even at 2 /im is still much 
smaller (^400Rs, where Rs = 2GMbii/c^) than at least the 
radius of the broad-line region^^ (^AOQORs), in the case of an 
untruncated multi-temperature black-body disk for our quasars. 
Significant geometrical effects will not occur unless the disk 
emission size becomes almost the same as the scattering region 
size. Therefore the near-IR polarized light spectra very likely 
reveal the intrinsic accretion disk spectra. 

The measured slopes, being as blue as the slope of the 
predicted shape F^, oc i^+^/^, strongly suggest that, at least in 
the outer near-IR emitting radii, the standard but hitherto unveri- 
fied picture of the disk being optically thick and locally heated is 
approximately correct. In this case, an implication is that other 
model problems at shorter wavelengths are associated with, or 
originate from, our lack of understanding of the inner regions of 
the same disks. We note that disk irradiation in limiting cases 
can contribute to the heating without changing the (y+^Z-^ spec- 
tral shape at long wavelengths; but it would not dominate the 
local internal heating in the outer radii considered here, and thus 
is not directly relevant except for some very specific cases^^-^**. 



The standard disk is also well known to be gravitation- 
ally unstable at large radii^^. These radii may well coiTespond 
to those emitting in the IR^" (^SOO-Rg for our quasars). In 
this case, if the disk is truncated at such a radius, the spec- 
trum will show a break, becoming even bluer at the longest 
wavelengths^^. Although statistically insignificant, our data do 
suggest that the near-IR slope is slightly bluer than the spec- 
tral shape Fy oc ly^^^^, with a hint of possibly becoming bluer 
at longer wavelengths. This can be followed up by extending 
the wavelength coverage with similar polarized Type-1 AGNs 
at lower redshifts. Such future measurements may pioneer the 
way to probe how and where the disk ends and how material is 
being supplied to the nucleus. 
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